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ABSTRACT 
A probabilistic model is developed for treating the 
problem of optimizing the intended trajectory and associated 
height of burst of a missile-warhead. The angle of re-entry 
is treated as a control variable. Computational techniques 


which may be operationally acceptable are described. 
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I INTRODUCTION 


Consider the problem of setting the fusing mechanism of a 
warhead aboard a ballistic missile in order to optimize the 
probability of destroying a target. 

The general problem as stated above becomes extremely 
complicated when one considers such things as missile in-flight 
failures, warhead failures, back-up fusing mechanisms employed, 
type of warhead and type of burst, distance from launcher to 
target (which will affect range and Azimuth errors), geodesic 
error, (i.e., to what degree of accuracy do we know the 
distance to target.), and more frustrating, the problem of 
partial kills, i.e., the partial destruction of the target or 
a portion of the target for a period of time less than the 
duration of hostilities. 

In order to reduce the problem to something that can be 
handled, it will be convenient to cover briefly three topics: 
the concept of ballistic coefficient, the reflected (mach) 
wave phenomenon and the definite range law (cookie cutter). 
Ballistic Coefficient: 

For our purposes, we may think of the ballistic coef- 
ficient as a parameter which describes how much a re-entry 
vehicle will be affected by the atmosphere (see Fig. 1). It 
is proportional to the ratio of the weight of the vehicle to 
its area, i.e., the square of the maximum diameter (see Ref. 4); 
e.g., for a given area, an increase in weight will tend to 
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"straighten'' the trajectory whereas a decrease in weight will 
result in a slower re-entry and thereby increase the time spent 
in the gravitational field. The result will be a more curved 
trajectory. 

Mach Effect: 

When a warhead is burst in the air close to ground (e.g., 
300 feet for a 1 KT warhead), there is a region of unusually 
high overpressures due to the merging of the incident and re- 
flected waves (see Ref. 1 and Fig 2). This region is known as 
the Mach Region. Because of this the curves in Fig. 2 are not 
monotonic but have the unusual shape exhibited there. 

Definite Range Law: 

This is the "cookie-cutter" concept that is frequently 
used. Because of its simplicity it seems to give useful 
results for many objectives. 

A couple of examples should be sufficient for our purposes. 
Consider a point target in space at which we are firing a pro- 
jectile with a kill range of one mile. Imagine the target to 
be at the center of a sphere of radius one mile. If we detonate 
this weapon anywhere within the sphere we consider the target 
destroyed; anywhere outside the sphere is considered a miss. Of 
course, no weapon behaves in such a manner. ‘''Partial kills" are 
generally the rule rather than the exception, and they are very 
difficult to quantify. The definite range law is a mathematical 
convenience. Presumably, this range can be regarded so that a 
controlled fraction of partial kills will be treated as misses 
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(and a controlled fraction will be scored as hits). 

As another example consider Fig. 5. If we imagine a target 
at (0,0) and estimate that it takes about 20 P.S.I1. of over- 
pressure for destruction, the 20 P.S.1I. curve shown could be 
considered as a "cooker-cutter"; i.e., a burst just inside 
will destroy the target and one just outside will not. 

In light of the above discussion, the following three 
simplifying assumptions will serve as the basis of this paper: 
1. The re-entry vehicle has a sufficiently high ballistic 
coefficient so that the trajectory may be assumed to be a 
straight line in the neighborhood of the target. Moreover, 
the vehicle is launched from a sufficiently great distance so 
that all possible paths may be assumed to be parallel in the 
neighborhood of the target. 

2. The overpressure contours of Ref. 1 (and hence of Fig. 2) 
are valid and the target damage is a function of overpressure 
alone. 

3. A definite range law holds and 20 P.S.1I. is necessary for 
destruction. 

The general objective is to select an optimum height of 
burst. Further examination of the problem reveals that some 
other parameters enter into the picture. The angle that the 
vehicle's trajectory makes with the surface can be selected 
in advance. There is an advantage to small angles since then 
the missiles are more difficult to detect (and hence defend 


against). Larger angles generally lead to higher kill 
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probabilities. The capability to vary this parameter will make 
defense more difficult. Given a trajectory angle, the optimum 
burst point becomes a function of height of burst and aim point 
(intersection of the intended trajectory with the surface). 
Locating the target at (0,0,0) and using the Cartesian system 
with 

\ = downrange distance 

X2= crossrange distance 

Nee vertical 
and introducing parameters 


~<J-= angle the trajectory makes with the surface 
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(see Fig. 3), the problem now may be stated: 


distance from aim point to target 


intended height of burst 


Given-> , choose the pair ( $ MA) so that the probability 
of target destruction is maximized. 

This paper must be viewed as a pilot study for solving the 
problem outlined above. It will be noted that a large amount of 
detailed work remains to be done in order to provide good answers. 
The methods herein can serve as outlines and guides. 

One of the more important results is the development of a 
graphical computation system using dividers, graphs of over- 
pressure functions and correction curves. Thus, the optimum 
solution may be obtained under operational conditions without 


the need of a digital computer. 





The organization of the report is as follows: The mathe- 
matical model and nature of the problems of mathematical 
analysis are presented in section Il. The computational problems 
are described in section III. The operational hand approxi- 
mations and corrective curves are given in section IV. The 
conclusions appear in section V. The Monte Carlo technique 
employed, along with the computer program, appear in Appendices 
A and B respectively. Appendix C contains a tabulation of 
the parameter values determined by the graphical method de- 


veloped in section IV. 
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CONSTRUCTION OF THE PROBABILISTIC MODEL; 
MATHEMATICAL CHARACTER OF THE PROBLEM 


We have already introduced the coordinate system; the 
angle the trajectory makes with the surface; the intended 
height of burst M ; the aim point ( oe. 0, 0). The trajectory 
is actually a random phenomenon and by assumption 1 can be 
characterized by the impact point and the incidence angle © . 
The actual height of burst is also subject to random errors. 
Thus, we let 

(Y1;, Yo, 0 ) be the random vector representing the impact 

point, 

(X,, Xo, X3 ) be the random vector representing the burst 

point 
A sample situation is given in Fig. 4. 

It is assumed that: 

1. Y1,¥9 has a circular normal distribution with mean vector 
a QO and covariance matrix cm where I is the identity 

matrix. 

2. %X3 has a normal distribution with sy oa variance g-“and is 

independent of Y1, and %. 


From Fig. 4 it is seen that X1 = Y, - X3 cot& 


2. Y2 


It follows that the probability density function of the 


detonation point is ' 
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The ballistic error parameter VU, , and the fusing error parameter 
O *should be determined from independent experiments. There are 
many types of warheads and fuses and we merely carry € , and (0; 
as "known" parameters. 

Since 20 P.S.1I. will destroy the target, the event of 
target destruction can be visualized in the (X1, X2, X3) space 
by rotating the curve in Fig. 2 about the vertical axis. If 
we call this set A, then the probability of target destruction 
may be expressed as the triple integral of the density function 
(2.1) over the set A 
aa P= JY (fen adr didy 

It may be wenn (eo clarify a small logical difficulty. The 
probabilistic model permits negative values for the height of 
burst. These are physically impossible and should be interpreted 
as a duded warhead that impacts the surface without exploding. 

The optimum solution can be obtained by taking the appro- 
priate partial derivatives of (2.2), equating them to zero and 


solving the resulting system of equations. Thus 
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Examination of the character of equations (2.4) is facili- 


tated if we make the change of variables 
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In this coordinate system the set A is transformed into a set B 
the exact nature of which is difficult to visualize (see Fig. 2), 
but we can say B is bounded since the transformation is linear. 
In the new system, equations (2.4) may be expressed as , 

ir Sf Ce- aly h t(4 4) | 

\ 27 [ 13 7 e = QO 
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Geometrically, the problem now may be viewed as follows: Given 


; { 





a set B, choose a centering point (5,0 ,M) so that the 
probability of B is maximized. 

All solutions of (2.6) are critical points which may be 
local maxima, minimum or saddle points. The nature of these 
points may be examined by means of second order derivatives. 
Letti Was V2 a A) 4 

etting oe - zL(Y,-S) tye(y-#) / 
\ = 
(2.7) A(t, %,%) = ra ¢g. 
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A critical point will be a local maximum if 


(2-411) “ays Xs 


It can generally be said that Daag (and oS , for the same 


2 
reason) is negative since the variance fhm"! AY yd ydhdy 
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of Y3 is one and (2.8) compares a truncated (by B) version of 


- Dae »O and Dam § 4 (Ret 3 pare ocy 


this variance with unity. The remaining part of (2.11) is a 
difficult question, however, and may require more specific in- 
formation about the set B. 

The set B plays a role in determining the number of critical 
points. For example, if B were a sufficiently long dumbell 
shaped region, there would be at least two solutions to (2.6). 
It would be interesting to examine the question of number of 
critical points if B were convex. With this condition it may be 
possible to show there is only one. This contingency was not 
examined since the set A (and hence B) is not convex. 

Since the set A is formed by rotation of Fig. 2 about the 
X3 axis, it can be shown that 
(2.12) Sf (x, £(X, 4, ¥,) dx, dx, dy, =O 
The same mn be true in the transformed coordinate system, and 
it follows from (2.6) and (2.12) that (€, 0fa) will be the 
center of gravity of the conditional probability density given 
CY]. .o a ommepelone to Batis F and (4 are unique solutions of 
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COMPUTATIONAL PROBLEMS 


Computer techniques are required for the solution of the 
system of equations (2.4), the examination of the critical 
points (2.11), and the determination of the optimum probability 
of target destruction (2.2). It was decided to use Monte Carlo 
techniques (see Appendix A) in this study, and the first step is 
to characterize the set A. The curve in Fig. 2 may be viewed 
as giving distance as a function. 

This curve is not given in analytic form. Also, note of 
height, i.e., g(X3 ) that there is a small interval of X, 
values for which g(X3 ) is double valued. It was decided to 
ignore the double valued feature and fit a polynomial of 
order 7 by the method of least squares, using 23 values read 
directly from the graph. This fit appears in Fig. 8. Using 
the polynomial as g(X3 ), the set A may be represented by 
5) Ee x1 4X2,K3 2 xi+ x5 < 5 k3)] 

The integrals in (2.2), (2.4) and (2.11) all are similar 
in nature. A rapidly converging iterative technique may 
generally be available for the solution of the system (2.3). 
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Let functions B1(M4, § ) and g2, § ) be defined so that the above 
equations are 
PA=Gi(M §) 5 FFF (MS) 

Then, using an initial approximation of M,,, a we can use 
the iterative scheme 
(353) Joa 7 Ga (Mn, i) oy, = 2, (Mon, ia 
provided it converges. 

The examination of the question of convergence usually 


goes along the following lines: If the series 
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numbers which will be the solutions of 2.4. This will take 


place if a constant r ( O€r<1) can be found such that 
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The. question a convergence of the scheme is more easily 
examined after making the transformation (2.5). Let this be 


done and write Es z IF OCn- 1 Pns dy, p)dy dy, I, 


and then 
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Proceeding as indicated, 
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From the theorem of the mean 
D(t)- O(r') < |a-r'| max PH) 


it is easily shown that 


/ / es 
(3.8) |g His = a < | 
(Wi S pee =F 
and this value may be used in 3.7. Thus we can write 


b (%, 4) - B%,-M,.,)| = P My -My-1| 

IBx- §,)- Bin 8.) PIF Fal 
Sharper bounds may be available if the maximum is 3.8 and the 
points ( r, ,O »M)? are constrained to be in the set B. The 
second term of 3.7 can pe estimated in a similar fashion: 

P| = ef, fGen 
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which can be treated as above. 
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Collecting, we have 
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which will have the form (3.5) if it can be shown that the 
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coefficients are less than one. This is a difficult analytic 
problem and is beyond the scope of this paper. The technique 
(3.3) was used in our computations and lend to answers that 
were stable in the light of our capability to compute by Monte 
Carlo. The signs of the determinants (2.11) were computed at 


the solution values and all were found to be positive. 
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IV 


A SIMPLIFIED APPROACH TO OBTAIN AN APPROXIMATE SOLUTION ; 
GRAPHICAL METHODS AND CORRECTION 


If the down range and cross range errors are assumed to 
be small, the problem can be reduced to a simple two-dimen- 
sional model since we intuitively feel that a ‘best’ trajectory 
is one that passes directly over target. 

At this point we will assume we follow our intended 
trajectory with probability 1 and find values of MM and 3 that 
optimize this conditional probability of destruction. 

If Fig. 2 were reflected about a vertical axis, Fig. 5 
will result. Again, assume the target to be at (0,0) and we 
wish to detonate our weapon inside the 20 P.S.I. envelope. 

We will make the following assumption to supplement those 
in section III: 

The warhead passes directly over target in a straight 
line trajectory making an incidence angle of with the surface. 
The collection of possible trajectories will intersect the 20 P.S.I. 
envelope as shown in Fig. 5. Denote the lower intercept(of a 
given trajectory) by a and the upper by &. 

The maximization takes place in two stages. First, holding 
© and 5 fixed, find the best value of AY as a function of eo 
second, vary y until the maximum is obtained. To implement this 


procedure, let 
- | ol ~on (Yr) 
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| ye: ao Je ie of Y 


13 





1000" 


\ 
x 
\ 
x 
x 
\ 
x 
~ 
. 
x 


MA 
fog 
Lt 





SACS SCPE 


Cah 
AN NSNSS 


ALLEANIAVALY ADTEVAD VE BUMOs Weavds f£ SALVA LIANES 


NAL 
‘ 
mi ow a ae Se 
8 *e 2 6 Ss. ic © «.. 
S 8 8 


9 
Ys 
~O 
CERT ILEE 


WHET NUE ETT 


ACCP SSCE 
ae | 
VCE SS 





WAY AISSS ACEH 
SOC ANNNe 1, 
| 

NYY ANNES 


ANAS 
AINAY CARY: | 


ae m x & 2 j 


lee 
‘DISTANCE FROM GROUND ZERO (FEET) 


oe 
| 


=< 


f 


3.67a. Beetiove erpr aba n the ground forget 1 “kilo ton bu rst eh pressu ure £. ang). 
; ; ans | laa ieee 


“Figur 


sieges SS = 


SAMPLE OF METHOD EMPLOYED IN ISOLATING OPTIMUM-INTENDED BURST POINT 





Taking the derivative and equating to zero 

= Short 
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Making the substitution t = Y-p 
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Since a function and its log have a miximum at the same points 


(Ke) = (447° 


Simplifying and solving for p 
q"-2emim = A-2amty” 
1 
Ve ee 


2(K-a) 
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The problem is degenerate for the case »& = a since the 
probability of detonation is zero. It is interesting to note 
that the above solution is independent of 0 and the velocity 
of re-entry. 


To verify we have a maximum, consider 
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Making the substitution t = Y - p and simplifying we obtain 
a 


after one integration by parts we see 
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which is true for all wx¥7Q 
The value of F at the point Mmax which is the probability 


of destruction given a particular line is given by 
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So far we have the optimum height of the IBP, but still 
must find out where to plan the detonation, i.e., to select €. 
Since we are interested in maximizing the probability of 
destruction (4.5), we see this is clearly equivalent to 
maximizing ( A -a) since the function is monotonic. Maximizing 
( X -a) is in turn equivalent to finding the longest line that 
can be drawn through Fig. 5 at an angle @. 

In Fig. 5, a series of parallel lines were drawn roughly 


every 20 feet at an angle@ to the horizontal and the longest 


ai 





line determined by the use of dividers, and choosing the one 

with max fo-a,f", If a maximum occurred between lines it 
ret 

was deemed to be at the mid-point of the two which gives an 

"accuracy' of 10 feet for a 1 KT weapon or 100 feet for a 1 MT 

weapon. 

By assuming the measurements were accurate in the sense 
that one can differentiate between a line being longer than an 
adjacent line (which is reasonable since the vast majority 
were easily discernable) we can say with probability one, the 

\/ 
maximum error is 10 W feet for a W KT weapon. 

Measurements were taken in this manner for several values 
of Gand the resulting values of fimax were plotted against © 
Readings were concentrated at points where the curve appeared 
to be rapidly changing. A similar procedure was used for values 
of @ vs X and (vs ia . in order to interpolate polynomials 


of degree 1, 2, ----- , 20 were fit in a least-squares sense to 


the data points. The degree chosen from the one with the 


smallest © [Y, - O, | f where Y: = emperical value 
Q 


and Oi = computed value. The above criterion was chosen since 
it was planned to interpolate values of G=5, 6, 7, ----- , 90 
and the measurements were felt to be fairly accurate to start 
with and large deviations in the fit would only tend to 
aggrevate the original data points. 

The polynomials were fed into the computer and evaluated 


at the points @-= 5, 6, ----- 90. The results were tabulated 





in table 1 where: 
jimax = planned burst point measured vertically along HOB axis 
XCOORD = planned burst point measured horizontally along 
"DISTANCE FROM GROUND ZERO" AXIS 
€ = Angle trajectory intersects ground measured from the 
horizontal 
ys Aim point measured positively to the right along "DISTANCE 


FROM GROUND ZERO" AXIS 


REMARKS CONCERNING THE GRAPHICAL APPROXIMATION 

The most interesting discovery is that the 'best' IBP is 
never over target; for low angles of re-entry, the burst is 
planned prior to target; for steep angles of re-entry, the IBP 
is past target. 

The sharp discontinuity in XCOORD at about 28° is caused by 
entry into the envelope in the region (850, 700) of Fig. 5 where 
the slope of the envelope and the slope of the trajectory are 
nearly equal. 

We will next apply the correction discussed in section IIl. 
The method developed above was used as the first approximation 
for the iterative technique developed in section III. The 


details follow. 


23 





| [n+(x,- Fx core) ]- Lo 


cry me F ee 
KEKE 9 (Ay) 


ff 
Pee (oh (ae F th TO) =0 


Re ft ala” es 


Ty 
m 
ct 
mre) 
WN 


=o 
T 


ee es [AE Ee [HF 


ENERN P+ Cot emn=e 
£,-MPz0 

Solving For AA and 5 
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let lal ) e, = First approximation obtained by using the 
method of section IV 


f E ; lo = BE, + Ma AG ieee 
/Am P ] M oe 

In general, about 10 iteratives were needed for@ > 10° 
and about 20 for@ = 10°. 

Data was collected for = 10°, 20°, ----- 80° for two 
sets of the standard deviation ( v v, ). Table 1 contains the 
percent of error in AA and f and Fig. contains a plot of vs 
percent of error. 


Thus, operational procedures can be developed that do not 


require a high speed computer. The first approximation is 
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d- = 500, cj = 200 6 = 200,71 = 500 


oa JA F Hh f 
10 t20 8 agis Jl +50 Zo 
20 Henri) +1 .34 +48 .9 tata 
30 + 4.28 -1.8 +127 Bee cae, 
40 ies 69 - .134 + See -24.2 
50 ee les 5 -4.8 - 6.15 -45 
60 see alOleys =2.0 =) 200 “41.5 
70 a 43 aS, = fan -65 
80 eo. So oe “oo “14.75 

TABLE 1 


PERCENT OF ERROR IN_“%/ AND ji BY USING FIRST APPROXIMATION 


yas 





JT = 200 
¢. = 500 
—O- M e PKILL 
10 453 2504 419 
20 431 1146 582 
30 432 719 652 
40 440 508 684 
50 438 369 692 
60 438 256 .699 
70 435 150 Ae 
80 438 84 .716 
g- = 500 
q, = 200 
& ye F PKILL 
10 365 2002 Oy, 
20 315 769 AO? 
30 391 658 .499 
40 459 661 569 
50 472 538 .615 
60 468 429 .637 
70 463 306 642 
80 471 155 649 
TABLE 2 


PROBABILITY OF DESTRUCTION AS A FUNCTION OF RE-ENTRY 
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obtained by the use of dividers applied to graphs like Fig. 5. 
The results are corrected by curves like Fig. 6. 

As was expected the method of section III gives a more 
accurate first approximation when the ballistic sigma is small. 
Probabilities of destruction can be predicted more accurately 
for higher values of €-. In general, the method of section IV 
gives values of Mand F that are low and nearer to target for 
small 4 (@ © 40°); high and further from target for higher ©. 

This has high intuitive appeal since if re-entry is at a 
small angle we tend to shoot for the point in Fig. 5 where the 
destructive envelope intersects the ground. But if ballistic 
errors are allowed, we run a high risk of having the vehicle 
impact with the earth before ever reaching target. Similarly, 
for steep angles of re-entry we tend to aim our errorless 
vehicle to pass thru the highest point on the envelope since 
the vertical distance traversed is a maximum. But if we allow 
errors, we run a high risk of passing to the left of the 
envelope of Fig. 5. 

Appendix C contains a list of starting points CA | fF? 
for @ =5, 6, 7, ----- , 90 along with associated destruction 
probabilities. If one chooses to use the curves given in 
Fig. 5, the destruction probabilities given may serve as a 
check on the final probabilities, the percent of difference 


decreasing as J, decreases. 
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Outside of Appendix C, all probabilities were computed 
using the polynomial fit to the curve since the curves were 
not believed to be accurate in the first place and we were 
interested primarily in technique. 

Graphs of the polynomial used and the actual curves are 


Found inerle. 72 
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FIGURE 7 


PLOT OF ACTUAL DESTRUCT CURVE AND FIT 


(both scales in hundreds of feet) 
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CONCLUSIONS 


The optimum probabilities of target destruction, and the 
optimum values of M and J appear as a function of © in 
Table 2 for parameter (9,6, ) =(200, 500) and (500, 200) ft. 
Notice that the fusing error parameter { is the more 
important parameter in terms of increasing the probabilities. 
This paper should be viewed as a pilot study. The method 
should be sharpened in the following ways: 
1. The overpressure envelope curves (e.g., Fig 2), should be 
established as a function of yield and mathematical methods of 
characterizing them (e.g., Fig 7), should be improved. 
2. Sharper numerical integration techniques should be applied. 
The present Monte Carlo method has a 3 standard deviation 
tolerance of + .015 for the probability of destruction. 
3. The technique should be applied to a large number of yields 
using realistic delivery errors and fusing parameters ( 9,andf ). 
4. The hand method for operational computations needs proper 
development, particularly in the area of sharpening the 
correction curves (e.g., Fig 6 and Table 1). 
For the parameter values used in this study, the probability 
of destruction was not very sensitive to changes in the values 
of M and Cee. This could change with other parameters. For 


example, a 5 MI weapon delivered from 2000 miles 


a3 





may well have a delivery error of say J, = 5 mites. If 50-2. scl, 
is required for target destruction, it may be necessary to 
determine AA and & very precisely. 

It was felt the assumption of a straight line trajectory 
is fairly realistic and that every moderately large deviation 
from a straight line will yield approximately the same results 
since the major dependency is simply in the fusing @ which is 
largely dependent on the particular warhead and not on the 
trajectory. Since curved trajectories mean lower speeds, the 
fusing error J will be smaller. It follows that the probability 
of destruction will be higher. Thus, the use of straight line 
trajectories will yield probabilities which may be used as 
lower bounds. 

The normalicy assumption in the fusing ¢ is felt to be 
pretty good if the warhead is fused by a timing mechanism or 
a radar, but barometric devices are notoriously non-normal. 

The circular normal distribution assumed for the trajectory 
is probably not realistic but is of little consequence since the 
model is broad enough to include different cross range and down 
range sigma's if necessary. 

The definite range law, as explained in the introduction, 
generally gives useful results and we feel it is quite 
appropriate for this problem. It may further be noted that 
this problem would not be solved were it not for this 


assumption. 
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APPENDIX A 


A brief description of the integration technique employed. 


Again consider a 20 P.S.I. curve in 3-space. We will first 
enclose the surface in a "box", then choose, at random, a 
coordinate inside this box. If the coordinate is inside the 
surface of revolution we record the "height" of the function 
at that point; otherwise, we record a zero. The sum of these 
observations divided by the number of observations will give 
the "average height" of the function over the surface. The 
average height multiplied by the volume of the box will 
approximately be equal to the integral of the function con- 
cerned. As the number of trials gets large the approximation 
will, of course, be better. 

If we consider each coordinate chosen as a Bernoulli 


2 
trial, the "worst" the variance could be is 0 = PS where 


N 
p = Py (coordinate is under surface) 
(as =P 
n = # of trials 


The standard deviation is Q 
In our case we chose n = 10,000 which yields an upper bound 


1 
(p = q = %) on the deviation ue ee 10,000 = 005 or a 


Be) acces! sof 1015 





We chose: 


Y1, Yo , Y3= Uniform (0,1) random variables 


Xj = 2 c(¥,- 0.5) x1 : Uniform €=8e, c ) 
Xo = 2. c(Y5- 0.5) X59 : Uniform (- c, c ) 


( b = 1000, c = 2000 in our case) 

If the triple of random numbers (X , X27, X3) drawn satisfy 

2 2 
the inequality X, == X,S g’ (KX) the function is evaluated at 
the point (X,], X9, X3), otherwise a zero is recorded. 
Let 
eg he ee 

h(X;) = integrand evaluated at (X1, X2, X3) if x) as (X3) 


— 10) 
we know that 


ee oe hn &) —> ELK (x)] = f{ Fe Por) AX 


VY — 00 
j= 


where P(X) = 1 in our case. The required integral is therefore 


approximated b a 
Le Z b oe n Sy “h (X;) 
1=0 


otherwise 


{- 


USE OF PROGRAM 
The following procedure is recommended for eeraeReaE fay, 210 PKILL: 
1. Function X GOF must be rewritten to fit destruction curves 
emp loyed. 
2. Values of § and@, must be changed in the main program to 


suit the weapon employed. 


oe) 
roy) 





3. Since only curves for a 1KT weapon are generally available, 
the problem must first be solved for this case using as a first 
approximation the value of Oars, found by the method of 


section IV and the parameters,4 § and@, must be multiplied by 
' U 


wo 


(¢ is not a function of the yield since it is inherent 


in the fusing mechanism). 


Conditional probability of destruction given the "best" 


trajectory when 


THETA = angle of re-entry 
A = lower intercept of trajectory with destruction envelope 
ALPHA = upper intercept of trajectory with destruction envelope 
Unease gevectical coordinate sous TbE 
XCOORD = down range horizontal coordinate of IBP 

7 = down range distance from target the trajectory 


intercepts the ground 
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THETA A XCOORD 3 ‘ 
5 199 999 *. 6595945 587 30.023 26550 
6 . £24,907 660.099 531 60.080 -44895 
iar : 389.185 653. 707 505 46.210 -~ 386 5 
BY 369. 360 645.386 49h (1.486) 3178 
9 51.793 610.607 487 47.129 ~2712 
10 JES). 999 636.991 479 77,2037 - 2382 
im 302. e523 055-2 OF Hos 50.5 =215n 
a 269.999 545 4 oh 92.987 -1923 
13 235.552 EN 304 438 91.381 ~1738 
i4 201.599 633.962 1420 Ot. 176 -1564 
AG 169.999 633.261 401956059 1399 
io tet. 37 YEA elim mae =124 ; 
7 114.999 631.038 369 111.540 -1096 
i8 89.067 630.126 356 alle tela. G64 
19 Si. 36 629.931 346 118.230 ~B19 
_20 29.999 630.912 218 Ani 20 -7 53 
ao””:*~=(<(t‘C OSS. ATS 56 Go ~ ter -6706 
22 | ~903000 6397948 333. 84.139 -620 
23 "00,000 644,412 334 - 75475 - 583 
2h, - 00,000 652.994 . 336 80.023 - 562 
he) 00.000 663.553 ee 2 WOE RS. - 5607 
20 | 153009) eS OL ee S00 
27 689.581 352 136. 230 -576 
28 704.321, Boo 139. 509 ~ 596 
29 719.641 366 =e. 566 -520 
30 7350102 0 37 -7% 982 ~ Sint 
onl 7506209 $00 Sms Tp Ss ve 
32 764. 844 388 -103.76&0 -687 
33 778.468 DDS =117.540 -704 
34 790.9657 ho2 ~131.339 _=-716 
£35 802.227 L408 ~i4h.830 | -725 
“30 Od 2.2207 ae ae “BES = 
‘37 i 0215092 24. -i70.220 -730 
* 38 828.780 427 -182.160 -728 
39 835.665 432 -193.820 ~7 24 
LO 84-1917 L377 ~205.430- , ~-719 
wa Ak: an a 2 
Lo § 53.488 WLS 7 229.370 -709 
43 859.259 1450 -241.960 ~-705 
Lb, 865.257 “454, 254.920 -702 
Ls 871.596 57 -268.120 -700 
~H6 873. 319 Loo =2016250 = =09 9 nnn Lc 
7] | | 885.401 462 ~293.940 ~700 
18 892.748 = 463 = 305.720 -700 
49 900.205 £65 ~ 316.090 ~700 
50 - | 907-564 466 324.550 -700 
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APPENDIX C (Con't) 


THETA A ALPHA nar XCOORD & 
oll 00.000 909.663 466 -330 695 
52 911.767 . 467 = 33 689 
59 | 921.521 » 467 = 33 681 
54 930.684 467 = 331 , 669 
OS 932.232 468 ~324 653 
56 933.000 468 -315 © 633 
Se 933.000 468 - 303 611 
58 467 -289 585 
59 : _ 467 -273 . $56 
60 467. -256 526 
61 am nl 467 =239 496 
62 466 =-223 ue 
63 465 -207 438 
64 146 -194 412 
65 L6 -183 389 
56 Loy -174 71 
67 464 -169 | 3 6 
68 464 -165 3845 
69 465 -164 | 338 
70 | 467 -164 . 333 
71 470 -164 "330 
72 a7 + -16 a 326 
73 479 -16 322 
74 48h -160 314 
75 490 -154 301 
7 94° =144 283 
77 495 -130 259 
78 . 49 -113 230 
79 ! 48 -937 196 
80 468 -7 32 160 
81 | LLO = 53 126 
82 401 -37 986" 
83 350 -269 801 
84 288 @~241 756 
85 i ed =297 865 
[om canna RH 00 C=) eee 
87 116 - 588 142 
88 124 -682 161 
89 221 - 538 136 
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: PROBABILITY OF DESTRUCTION 
SIGMA = 10 50 100 200 500 1000 


1.009 1.000 1.000 .978 .640 °.353 
6979 643 355 

A 0979 646 5S ole 
~ 980 648 03538 
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